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Abstract. We compute a variety of operator-operator correlation functions to third order in the 
MS scheme in the chiral limit. These include combinations of quark bilinear currents with gauge 
invariant operators such as moments n = 2 and 3 of the flavour non-singlet Wilson operators of 
deep inelastic scattering and moment n = 2 of the transversity operator, as well as the correlation 
functions of the latter operators with themselves. The explicit values of these gauge independent 
correlation functions are required to assist with non-perturbative matching to lattice regularized 
calculations of the same quantities. As part of the computation we determine the mixing matrix 
of renormalization constants of these non-singlet currents with their associated total derivative 
operators at the same twist to three loops in MS. Such operators are crucial in extracting 
renormalization constants for the operator-operator correlation functions which are consistent 
with the corresponding renormalization group equation. As a by-product we deduce the i?-ratio 
for the tensor current to third order in the MS scheme. 
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1 Introduction. 



Lattice regularization of non-abelian gauge theories has provided many insights into the non- 
perturbative properties of the quantum field theory underlying the strong force, known as Quan- 
tum Chromodynamics (QCD). For instance, now that dynamical quarks can be treated via the 
use of powerful supercomputers, the meson and hadron spectra are on the whole in remarkable 
agreement with experimental data. One current problem of interest for lattice regularization 
is the measurement of Green's functions relevant for, say, deep inelastic scatttering. For in- 
stance, whilst the perturbative renormalization of the underlying twist-2 flavour non-singlet and 
singlet operators are known to three loop accuracy in the MS scheme, [H [21 [3l HI O [6l [7] , the 
associated matrix elements are also required but they can only be fully computed using non- 
perturbative methods. These are crucial to fully understanding the structure functions of the 
original hadrons and mesons which are broken up in deep inelastic experiments. Such matrix 
elements are, however, accessible via lattice regularization with notable progress through the 
years via collaborations such as QCDSF, [H [9l [TOl dH 111 13] , and others, [H [ISl H HZ] . How- 
ever, such measurements of Green's functions on the lattice need a variety of techniques to allow 
comparison with continuum results. Therefore, we will first briefiy discuss some of the relevant 
issues in a general context before formulating the aims for the current article. 

The first issue is that the lattice computations necessarily renormalize their operators and 
Green's functions using a renormalization scheme which is not the standard MS one. The general 
name for the scheme we refer to in this context is regularization invariant (RI), [IB]- Though in 
practice the main scheme used is referred to as RI'. (In some articles this is synonymous with 
RI-MOM.) Therefore, one needs a means to convert lattice results from RI type renormalization 
schemes to the reference scheme of MS, [18\ I19| . For certain classes of Green's functions the 
continuum definition and use of RI' (and RI) have been given in three and four loop renor- 
malization of (massless) QCD in arbitrary covariant gauges, [20 \ I21[ [22| I23j. The second main 
issue to deal with is that of matching lattice regularized results for the Green's functions with 
the corresponding continuum results. There are two main approaches for this. The first is the 
Schrodinger functional method, [Ml EH], upon which we make no further comment. The second 
is the matching to explicit perturbative QCD results for the same Green's function which is 
what we concentrate on here. In this approach, [18\ [T9l [20l [2T1 [22l I23j . the main ethos is to 
compute the Green's function to as high a loop order as possible in perturbation theory as a 
function of the gauge coupling constant, g. Then the lattice computation should match on to 
the continuum behaviour in the same renormalization scheme in the high energy limit. Having 
a perturbative result to as high a loop order as is calculationally feasible will in principle lead to 
a more accurate extraction of numerical values for the Green's function in the non-perturbative 
range of interest for, say, nucleon structure functions. 

As already indicated, previous work in the continuum concentrated on a variety of gauge 
invariant twist-2 fiavour non-singlet operators, O, and the perturbative evaluation of the flavour 
non-singlet Green's function {iJj{p)O{0)^{—p)) in the chiral limit, [20^ [2T1 [22l I23j . where p is 
the momentum flowing through the Green's function and ip is the quark field. The operators 
considered for this Green's function were the Wilson and the transversity operators to and 
including moment n = 3 and 4 respectively and various quark bilinear operators such as the 
tensor current, [211 122^ [23]. However, whilst such results were useful in many ways, they suffer 
from one major drawback. This is simply stated by noting that although each of the operators 
considered was gauge invariant, the Green's function itself was gauge dependent. Although 
ultimately one was only interested in the Landau gauge, the results of [HI [22l [23] were provided 
in an arbitrary linear covariant gauge. Whilst this was not too problematic for continuum 
calculations, from the point of view of lattice regularization one has also to fix the Landau 
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gauge. However, this is a tough exercise in itself and in principle could open up reliability 
issues to do with say ensuring the Gribov problem was avoided. To circumvent these lattice 
regularization gauge fixing issues another approach has been devisecfl- Briefly to extract the 
appropriate renormalization constants for the operators and hence determine the finite parts 
of the Green's functions, the approach is to consider gauge independent correlation functions 
of gauge invariant operators. In this way the potential gauge fixing ambiguity never becomes 
an issue in the first place since the gauge does not then need to be fixed on the lattice. More 
specifically the appropriate correlation function to consider is {0{p)0{—p)) . However, for, say, 
high moment Wilson operators the increase in the number of covariant derivatives may lead to 
too noisy a numerical signal for extraction of meaningful values of the Green's function. Hence, 
rather than consider this diagonal correlation function, a simple proposal would be to analyse 
an off-diagonal correlator such as {O^ {p)0'^ {—p)) where 0^{p) is a Wilson operator, say, and 
0'^{p) is a simple quark bilinear current operator. This has to be chosen in such a way that the 
Green's function is not simply trivial in the chiral limit. However, information on the Wilson 
operator renormalization constant can still be derived. 

Having reviewed the background and key issues we now indicate the main aim of this article. 
It is to simply to provide the explicit values of the appropriate operator correlation functions, 
{O^ {p)0'^ {—p)) , relevant to the lattice problem to as many loop orders in perturbation theory 
that is calculationally feasible. Specifically we will focus on three sets of flavour non-singlet 
operators used in deep inelastic scattering. These are the Wilson operators with moments 2 
and 3 and moment 2 of the set of transversity operators. Several correlations of the operator 
with itself will be provided as well as the appropriate non-zero off-diagonal one. We will present 
results to third order or three loops in the MS scheme as a function of the momentum flowing 
through the 2-point function. We note that whilst we compute three loop Feynman diagrams, 
since it is clear that the leading diagram of the correlator is independent of the strong coupling 
constant, a, then the results will be to 0{a?) inclusive where a = / {IQ-k'^) . However, it will 
also become evident that it is not possible to consider correlators of higher moment operators 
and expect to evaluate the correlation functions to the same three loop order. Whilst our main 
motivation is to provide the flnite parts of these correlation functions several technical issues 
need to be addressed to obtain the correct answers. For instance, there is an assumption that 
the flavour non-singlet Wilson and transversity operators do not mix under renormalization. 
It will turn out that this observation needs to be clarified within the present context. Their 
three loop MS anomalous dimensions are available, [H El |3l HI El El [7], and the lower loop 
results were originally obtained by considering the renormalization of ('i/'(p)O(0)V'(— p))- In this 
momentum configuration the mixing is not relevant. However, in the momentum configuration 
for the correlators of the present article {0^{p)0'^{—p)), since a momentum fiows through the 
operator the mixing is relevant and cannot be neglected. Suffice to say at this point that the 
additional operators are gauge invariant but total derivatives. Therefore, as part of our correlator 
renormalization programme we have had to compute the relevant anomalous dimension mixing 
matrices to allow us to extract consistently renormalized correlation functions. Such results will 
no doubt be important for other areas of deep inelastic scattering such as generalized parton 
distribution function analyses. 

The article is organised as follows. Section 2 introduces the notation, operators and general 
features of the correlation functions we consider throughout. The general renormalization prop- 
erties of the underlying correlation functions are discussed in section 3 including the operator 
mixing issue. Section 4 is devoted to the very mundane but important exercise of recording all 
the results for the Green's functions we have considered. As a spin-off we record the i?-ratio 
for the tensor current to third order in section 5. Finally, after concluding remarks in section 6, 

*The author is grateful to Dr P.E.L. Rakow and Dr R. Horsley for their patient enhghtment on this point. 
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several appendices are provided. The first records the Lorentz tensor decomposition of several 
operator correlators. This is necessary since the lattice requires the use of operators with un- 
contracted indices. In the renormalization of the matrix elements {ip{p)O{0)'ip{—p)) of [HEIIS], 
the Lorentz indices were contracted with a null vector with = 0. This was because, for 
example, the Wilson operators were traceless and symmetric and this contraction excluded the 
part with metric tensors in order to ease the extraction of renormalization constants directly. 
Here, since we will use a very specific computer algebra package and algorithm which only op- 
erates on scalar Feynman integrals, we need to project out the relevant scalar amplitudes with 
respect to some tensor basis, which is, of course, not unique. This is discussed in Appendix A. 
The remaining appendix records the explicit numerical values of the various finite parts of the 
correlation functions for the colour group SU{3) which were originally presented in exact form 
to O(a^) in section 4. 



2 Preliminaries. 

In this section we define our notation and operators and discuss the operator correlation functions 
from a general perspective. Throughout we use the standard QCD Lagrangian with massless 
quarks to immediately put us in the chiral case and an arbitrary linear covariant gauge fix- 
ing which is parametrized by the (renormalized) parameter a. However, since our correlation 
functions involve gauge invariant operators and therefore are gauge independent, a will never 
actually appear in any of our final correlator expressions. Though we stress that at no stage have 
we set a = internally in our computations. Its natural cancellation is a strong internal con- 
sistency check on the construction of, say, our Feynman rules and the operator renormalization. 
Given this we define the general correlation function as, [26 \ 127 1 [28], 

= J d'^^e^^^(0|O;.,...,„Xx)O^.,....„X0)|0) (2.1) 

where q is the momentum (with q'^ = — Q"^) and we have labelled the Lorentz indices of the 
respective constituent operators by a different parent Greek letter for clarity. At this point it is 
worth noting that we closely follow the procedures of |261 [271 128] which are excellent reviews of 
calculating ()2.ip for quark current operators. The Green's function itself is illustrated schemat- 
ically in Figure 1 with the momentum flow made explicit. Our notation needs explanation. We 
use superscripts i and j to denote the left and right operators O of the correlation function as 
indicated in Figure 1 where the momentum flows into the left operator and out through the 
right operator. In principle, these operators are different whence the two distinct sets of Lorentz 
indices {/ii} and {vj}- Though for some cases we will include the quark mass operator which 
has no Lorentz tensor structure and so no such indices will be formally required. The flavour 
indices have been omitted to avoid cluttering the notation further. It is understood that there 
is a flavour generator included within each operator and later we will note the internal account- 
ing method used to ensure that we derive results only for the correlation of flavour non-singlet 
currents as opposed to flavour singlet currents. For the latter there would be an additional but 
different operator mixing problem from that into total derivative operators. This singlet mixing 
is already well documented, [Il[2l[3]. In other words flavour singlet quark blinear current oper- 
ators can mix into purely gluonic operators with the same twist and quantum numbers, as well 
as gauge variant operators with the same properties but constructed from, say, Faddeev-Popov 
ghost fields. Moreover, one would also have to handle equation of motion operators too. We 
mention this aspect for completeness and note that as far as we are aware there is currently no 
lattice proposal to examine the flavour singlet case and we therefore will not devote any time to 
it here in the analogous continuum problem. 
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Figure 1: Operator correlation function {0^{q)0^ {—q)). 



Since we will consider a variety of operators we introduce a shorthand notation, akin to 
[26\ [271 128j . for the superscripts i and j to indicate which operator appears in (j2.ip . These are 
listed below as 



s 


= Tplp 


V 




T 




W2 




dW2 


= Sd'^i^Yi') 






dWs 




ddWs 


= 59^5'^ (^7'"V) 


T2 




dT2 





where S denotes the appropriate symmetrization in the Lorentz indices as well as the operator's 
tracelessness which we will discuss shortly. This list already reveals our hand in terms of the 
operator mixing issue. Moreover, in choosing this notation we have disguised to a degree what 
some of the total derivative operators are. However, we choose to work within certain sectors. By 
this we mean, using the Wilson operator of moment 3 for illustration, that W3 is the main label 
for that sector as well as being the parent operator. In its renormalization it spawns offspring 
total derivative operators denoted by one or more d symbols. These are labelled dW^ and ddW^. 
However, clearly from (j2.2p these derive from the parent W2 operator of the preceding or lower 
sector and the vector current of the first sector. It therefore might have been apt to choose 
the respective notation for these to be dW2 and ddV for labels. We have chosen not to do so 
because for the W3 sector the number of Lorentz indices are the same on each combination of 
these operators. Moreover, one would have to count the derivative labels to deduce which sector 
the label was associated with. Hence the Lorentz projection tensors used to project out the 
scalar amplitudes we will compute are the same for these combinations. So for the W3 sector 
there is only one basic projector. Whilst this means we will use W3 as both an operator label 
and sector indicator, this we believe will minimize the confusion in choice of notation when 
trying to ascertain which sector, say, dV belongs to especially when transversity is dealt with 
in parallel calculations. 

Next in (12. 2p we choose the symmetrization with respect to Lorentz indices and the traceless- 
ness conditions in the standard way. For transversity, which involves the 7-matrix commutator 
a'^'^ = [7^,7^^], the definition is not the same as for the Wilson operators, \29 \ \30 \ [31]. For the 
three moments we consider, the explicit definitions of the symmetric traceless operators are, in 
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d-dimensions, 



S fiua 



where 



1 

fiua 



1 
6 



{d + 2) L '^'^ •S' '^P 'lua'^s Atp 'IcT^^s vp 

2 



(2.3) 



(2.4) 



which are clearly consistent with the definitions in [2H [22| I23j . We note that 05'^,yo- is the 
intermediate definition of the symmetrized operator. As in [21^ [22l [23] we have derived the 
d-dimensional versions. Although the lattice computations are in strictly four dimensions, we 
will dimensionally regularize in d = 4 — 2e dimensions where e plays the role of the regularizing 
parameter. The renormalization constants will have a Laurent series in e when subtracted in 
the MS scheme. The main calculational tool is the use of the Mincer algorithm derived in [32]. 
It is ideal for the present work as it evaluates massless 2-point Feynman diagrams to the finite 
part at three loops in dimensional regularization. The correlation function ()2.ip clearly falls into 
this category and does not require infrared rearrangement or external momentum nullification. 
More practically the Mincer algorithm has been encoded in the powerful symbolic manipulation 
language Form, [33j, in [M], which therefore allows for a fully automated computation. For 
instance, once the Feynman rules for the operators are derived consistently then the algorithm 
is applied to produce the finite parts. Crucial to this is the electronic generation of the Feynman 
diagrams via the Qgraf package, [35]. These are then converted into Form input notation for 
application of the MiNCER algorithm by systematically including the Lorentz and colour indices 
for the gluon, quark and Faddeev-Popov ghost fields. For the present calculation there are 1 
one loop, 8 two loop and 109 three loop Feynman graphs to be evaluated in principle for every 
combination of operators in ()2.2p we consider here. These totals include graphs where there are 
two covariant derivatives in each operator as occurs for the W3 sector. In the Qgraf generation 
of graphs we restrict the diagrams to the one particle irreducible, no tadpole and no snail set-up 
since we are dealing with massless fields. So, for example, there are no closed gluon loops at the 
location of an operator. Such graphs are trivially zero in dimensional regularization but may arise 
in, say, a lattice regularization. That aside, when, for example, ddW^ is part of the correlator the 
majority of the 109 three loop graphs will in fact be trivially absent. For practical purposes it is 
best to have the most general set of diagrams for the full calculation rather than design Qgraf 
routines for specific cases and potentially omit graphs which contribute. By the same token 
the presence of the covariant derivatives in the VF3 sector means that the explicit evaluation is 
slowed significantly on the available computers. Therefore, when this occurred we chose to run 
each Lorentz projection individually in series, which improved run times substantially. Given 
this we note the final aspect of our notation and that is the decomposition of the correlation 
function into the explicit scalar amplitudes, n*^-^j(g). These are defined by 



'-^lil...p.„.ui. 



k=l 



(2.5) 
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where ij. v -}(^) ^'^^ Lorentz projectors. The subscript (k) (and also (/) later) 

label the sector. How these projectors are derived and their explicit forms are relegated to 
Appendix A. However, we note the number of projectors for each of the eight correlation function 
sectors we focus on here, njj, is given in Table 1. Clearly the number of projectors increases 
with the number of free Lorentz indices. Given this decomposition then to find each individual 
scalar amplitude n|-^^(g) we multiply (j2.5p by the appropriate element of the inverse projection 
tensor which is defined for each sector in Appendix A. It is then this Lorentz scalar object which 
is put through the MiNCER algorithm. 





S,S 


v,v 


T,T 


V,W2 




W2,W2 


W3,W3 


T,T2 




1 


2 


2 


2 


2 


3 


4 


4 



Table 1. Number of projectors for each correlation function sector. 



As is usual with a renormalizable quantum field theory each amplitude is divergent. To 
extract the explicit divergence we follow the algorithm of [36] derived from automatic multi- 
loop computations. In general terms, one computes the Green's function as a function of bare 
parameters. Then the renormalized variables are introduced by the simple rescaling definition. 
For example, for the bare and renormalized coupling constants go and g respectively, we use 
go = gZg where Zg is the coupling constant renormalization constant. Though in dimensional 
regularization we will use go = fJ-'^gZg where /i is the arbitrary renormalization scale present due 
to the regularization. Here the explicit renormalization constants for the operator correlation 
functions are derived and discussed at length in the next section where the associated renormal- 
ization constants are also constructed. It suffices to say at this point that there is an extension 
to the algorithm of [36] in that the bare operators in the Green's function have to be rescaled to 
their renormalized operator without neglecting the mixing into other operators. Moreover, the 
explicit forms of the gauge independent operator correlation function renormalization constants 
are equally as important as the finite parts of the correlators from the point of view of assisting 
with the matching of lattice results to the corresponding continuum values in the high energy 
limit. 

There is one concern with relation to the topology of the graphs which needs to be addressed. 
That is that we need to ensure that within our computer programmes we are in fact calculating 
the correlation functions of flavour non-singlet operators as opposed to flavour singlet ones. For 
instance, without the presence of a flavour matrix of some sort graphs which have a closed quark 
loop and only include one of or but not both together must be set to zero. If not it would 
be a contribution to a flavour singlet operator correlator. Therefore, whilst we have formally 
omitted flavour indices in the definition ()2.ip . our Feynman rules for the operators actually 
include a flavour matrix for each operator. Denoting this by A* where i labels the left operator 
1 or right operator 2, then at an appropriate point of the computation terms with tr (A^) tr (A^) 
are set to zero to only leave terms proportional to tr (A^A^). This is then formally set to unity 
since it flags the flavour non-singlet contribution uniquely. 

Finally, we comment on how we have chosen the correlation functions presented here. The 
set we consider is listed in Table 1. First, our choice is motivated by ensuring that the corre- 
sponding lattice calculation has a minimal set of covariant derivatives to handle. Second, we are 
constrained by the masslessness of the problem. For instance, as is clear from Figure 1 and (j2.ip 
we are dealing with closed quark loops. Therefore, one must have an even number of 7-matrices. 
As the quarks are massless and each quark propagator has exactly one 7-matrix the sum of 
7-matrices present in both operators of (j2.ip has to be even. Therefore, whilst from a lattice 
point of view it would be simple to have the quark mass operator as the off-diagonal element for 



7 



the W2 and VF3 sectors this correlator is trivially zero. In other words in the presence of quarks 
with generic mass niq then the correlator will vanish as 0(m|) where > 0. Hence, for W2 and 
W3 they have to be paired with V. Similarly as T2 involves a^'^ one requires an even number of 
7-matrices for the other operator of the correlator. Naively one would assume that this natural 
pairing would be with S. However, that leaves free Lorentz indices only on one operator and 
for T2, given the symmetry properties of the operator itself via a^'^ , it is not in fact possible 
to decompose the correlator into Lorentz tensors built from the metric tensor, ij^'^ , and the 
momentum, q^. Therefore, we have had to pair T2 with T. Whilst these off-diagonal operators 
will probably be the ones of most interest to lattice computations we have chosen to consider 
the diagonal sectors {W2,W2} and {VF3, W3} as well. There are several reasons for this. With 
a second avenue available to extract information on all the W2 and W3 sector renormalization 
constants used for the lattice, these will actually give useful consistency checks provided suffi- 
cient computation power is available for the lattice calculations. Next, the operator mixing issue 
is a novel feature of these correlators and we choose to consider them to ensure that we have 
obtained the correct overall picture of view from a calculational and renormalization point of 
view. A final, less firm, motivation is that the diagonal correlators of quark bilinear currents are 
useful to derive decay rates via the i?-ratio formalism. (See, for example, \26 \ \27\ [28].) Whilst 
those for W2 and W3 are tenuous in this respect, and we take them no further than finding the 
amplitudes, we do evaluate {T, T} for this reason and construct the corresponding i?-ratio for 
the tensor current as a by-product of our full computation. 



3 Renormalization group. 

In this section we concentrate on general aspects of the renormalization of the operators we are 
interested in, ()2.2p . and the construction of the renormalization group equations satisfied by the 
renormalized operator correlation functions. There are essentially two parts to this. The first 
relates to the operator mixing which is a separate exercise unrelated to the correlation functions, 
whilst the second is an application of the mixing property. 

For the operator mixing issue we note first that the quark current operators S, V and T are 
clearly mutiplicatively renormalizable in the chiral limit and mass independent renormalization 
schemes. Therefore, we concentrate on the three sectors W2, W3 and T2 and justify our choice 
of operator basis. The first and last of these are similar, so we will specifically consider W2 
and W3. First, for W2 we note that the usual Lorentz symmetric and traceless twist-2 flavour 
non-singlet operator used in deep inelastic scattering is 

= Sijj^D.i; . (3.1) 

This is not independent since one can of course add the independent operator S (D^ip) 'jf^'ip to 
the set of operators with the same symmetry properties. However, we want to make use of known 
renor mahzation results for O^J', [HElElllElElElETllaH], and using the latter noncanonical 
operator, whilst not difficult from a technical point of view, is not the only operator independent 
of it. Instead the operator 

O;^^ = Sd^ (V^T.^) (3.2) 

is independent of W2 and with dW2 and W2 we can obtain S (D^ip) j^ip as a linear combination. 
Likewise for the T2 sector the analogous basis is 

0%l = Sd^{i;a,^i;) . (3.3) 
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At the next sector higher one now has an object with three Lorentz indices and so one would 
expect three independent operators. Again we wish to retain the renormahzation properties of 
the parent operator which means choosing the basis as {W3, 5W3, SSVFs} where 

Of^^ = Sd.d^iijj^i;) . (3.4) 

Clearly the latter two are total derivatives of the set {1^2)91^2} and if one were to extend 
to the next sector level that set would involve the parent and the total derivatives of the 
VF3 sector. From the explicit calculation of the operator anomalous dimensions there is an 
important computational advantage from choosing the basis in this way. Though it should be 
stressed that at higher levels the choice of basis is arbitrary and one could in principle choose, 
say, Sdfj, [(^Dy'il)) "faip) as an independent member of the set. 

With the choices we have detailed for each sector, there is mixing under renormahzation but 
clearly each mixing matrix of renormahzation constants, Zfj, is upper triangular where 

Ooi = ZfjOj (3.5) 

relates bare operators, denoted by the subscript o, to their renormalized version. To be explicit 
the matrix for the W2 and T2 sectors is 

'S - if I ^ p-«> 



and that for W-^ is 



yO yO yO 
^11 ^12 ^13 



= \ Zg Zg I . (3.7) 
V Zg 

Here we have chosen to simplify our notation by using numbers to denote the mixing matrix 
elements rather than the more clumsy {W2, {W3, SWs, ddW^} or {T2, 8X2} as subscripts. 

Given these matrices we then define our anomalous dimension mixing matrix elements, ")gj{cL), 
formally as 

7? = /^l^ln^^ (3.8) 

where 

Here P{a) is the /3-function and 7^(0, a) is the anomalous dimension of the gauge parameter 
where we follow the conventions used in [21] to define its renormahzation. Although all our 
renormahzation constants will in fact be independent of a we have included it in (jS.Op as it 
technically appears as a formal parameter in the QCD Lagrangian. For a renormahzation where 
there is operator mixing, (j3.8p is invariably given as the formal definition of the anomalous 
dimensions. However, for practical purposes in the derivation of the operator correlation function 
anomalous dimensions it is more appropriate to give the explicit consequences of p.Sp . For 
sectors 14^2 and T2 we have 

d 



= 7?i(a)^fi + /x^Zg 



= 7n(a)Zf2 + l?2{a)Zg + fi—Zg 



d 



l?2{a)Zg + fi—Zg (3.10) 
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and for we have similar relat: 



ions, 



= 

= 

= 

= 

= 

= 



7S( 

y?i( 

.Pi 



^11 


+ 






yU 


+ 


7f2(«)^2''2 




yLJ 
^13 


+ 


7f2(«)^2''3 


+ 7?3(a)^3''3 


yO 

^22 


+ 


^d^^22 




yO 
^23 


+ 


72''3(«)^3''3 




yO 
^33 


+ 







(3.11) 



t P1+P2 




Figure 2: Green's function, pi — P2)'^{P2)), used to renormalize the operators O*. 

We now turn to the practical problem of evaluating the anomalous dimensions explicitly to 
the loop order necessary for the operator correlation function renormalization at three loops. 
Whilst this is actually 0{a^) we have determined the mixing matrices to O(a^) partly for com- 
pleteness and checking reasons but also because of their potential use in phenomenological prob- 
lems where a momentum flows out through the operator itself. For the mixing matrix anomalous 
dimensions the main issue is the determination of the off-diagonal elements. In the original ap- 
proach of [U [21 [3] the operators were inserted in the Green's function (■0(p)O(O)V'(— p)) whose 
more general version is illustrated graphically in Figure 2. As the operator is a zero momen- 
tum insertion the contribution to the renormalization of the off-diagonal part of Zfj cannot 
be deduced as the total derivative operator insertions vanish trivially for this momentum con- 
figuration. So only all the diagonal elements of our mixing matrix can be deduced via this 
momentum routing. There are two remaining choices for routing momenta with one nullifica- 
tion, which is necessary if we wish to apply Mincer as our tool of computation. With the 
choice {ip{p)O{—p)ip{0)) it is clear that one can access the elements Z^'^, ^'^'^ ^23^ once 
the respective values for Z^'^, Zjl and Z^^ have been included from the original results of 
[HIllEliaElElEllMlllQlil]. The determination of Z}^^ and Zj^^ is more difficult. This is be- 
cause there is only one momentum choice which makes contact with both these renormalization 
constants at the same time. We circumvented this difficulty to two loops by not omitting the 
terms involving ln(p^//i^) which derive from the dimensionality of the loop integrals at each loop 
order. One ordinarily ignores such terms in a renormalization since they cancel trivially in a 
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renormalizable theory. Retaining them means that divergences such as ^ ln(p^//x^), ^ ln{p'^ / fj?) 
and ^ (ln(p^///^))^ are present at various loop orders but some of their coefficients involve coun- 



terterm parts of Z^2^ and Z^g^ in their coupling constant series and Laurent expansion in e. 
They give extra constraints on zj^^ and which allowed us to decipher the sector mixing 
matrix to two loops. At three loops we can only derive a simple relation between the simple 
poles of each renormalization constant. However, this is not in fact necessary for our operator 
correlator results. Again the main tool is the Mincer package written in Form, |35l [36l I37j . 
The graphs are again generated by Qgraf, [3H]. For the maximum number of possible covariant 
derivative terms that can arise, there are 3 one loop, 37 two loop and 684 three loop Feynman 
graphs to evaluate. As before these are one particle irreducible graphs without snails or tadpoles. 

Given these considerations it is evident that sector W4 would be more difficult to deduce 
completely to three loops due to a similar momentum routing issue. However, it could be 
determined to two loops in principle by performing a four loop calculation with ln(p'^//u^) terms 
included. If there were a four loop Mincer routine available then this would be a viable 
proposition. Therefore, in principle, the procedure to determine the mixing matrix for the non- 
singlet sectors Wn is available but in practice requires the computational machinery to evaluate 
the off-diagonal mixing matrix elements. All that remains is to record the explicit values which 
are all given in the MS scheme. First, for completeness and for comparing with the conventions 
of previous calculations the vector and tensor current anomalous dimensions, \39 \ \iO \ are 



7^ (a) 



Cpa + [257Ca 



UICf - 52TFNf] 



18 



+ 



13639C7| - 4320C(3)Ci + 12096C(3)CaCf 

- 20469CaCf - 1728C(3)CArpiV/ - mQCATpNf 

- 6912C(3)C|. + 6570C72 + 1728C(3)CFTpiV/ 

Cpa^ 



108 



(3.12) 



where ({z) is the Riemann zeta function, Ca and Cp are the usual colour group Casimirs and 
Tp is defined by Tr {t"'T^^ = Tp5°'^ where are the generators of the colour group. For the 
two sectors with two operators we have 



7i?(a) 



Cpa + 

1 
243 



1 

27 



376CaCf 



112C^ 



(5184C(3) + 20920) C\Cp 



USCpTpNf c? 
(15552C(3) 8528) CaC\ 



- (10368C(3) + 6256) CACpTpNf + (10368C(3) - 560) C|, 
+ (10368C(3) - 6824) C^TpNf - 896CpT^N^] + 0(a^ 



71^^ (a) 



56C| 



+ 



-Cpa + — 
3 27 

^ ^(7776C(3) + 4264) CaC| 



243 



+ 



188CaCp + GACpTpNf 

(2592C(3) + 10460) C^Cp 
(5184C(3) + 3128) CaCpTpNj - (5184C(3) 



(5184C(3) 



3412) ClTpNf + AASCpTjN^ 



280) C| 



72^^^ (a) 



(3.13) 
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and 



SCpa + - 35CaCf - 9Cl - UCpTpN^ 



+ 



108 



12553C1Cf - 7479CaC|, + 1782C|. - (5184C(3) + 4168) CACFTpNf 
+ (5184C(3) - 3240) C^TpNf - SGSCpT^N^] + O(a^) 



j'(^(^a) = - Cpa + hsCpTrNf - 45C|^ - 29CaCf 

18 L 



+ 



108 



18 

(6048C(3) - 6495) CaC^ - (2160C(3) - 543) ClCp 
+ (1728C(3) + 76) CaCfTfNj - (3456C(3) - 2394) 0% 

- (1728C(3) - 2208) C^TpNf + lUCpT^N^] + O(a^) 

Cpa + ^ [257CaCf - lllCl - ^2CFTFNf\ 

+ t4t [(13639 - 4320C(3)) C\Cf + (12096C(3) - 20469) CaCI 

- (1728C(3) + 4016) CaCfTfNj - (6912C(3) - 6570) 
+ (1728C(3) + 1176) CJ^TfA^/ - 144CFT|Ar/] + O(a^) 



108 



(3.14) 



Given the issues with computing the full set of anomalous dimensions to three loops for the W3 
sector we record that the results are 

7n'(«) = y C^a + ^ [8560CaCf - 2035C|. - 3320CFrFAr^] 

(285120C(3) + 1778866) C^Cf - (855360C(3) + 311213) CaC], 
- (1036800C(3) + 497992) CaCfTpN^ + (570240C(3) - 244505) Cl 
+ (1036800C(3) - 814508) C^TFNf - 82208CFTj,N^] + 0{a^) 

-Cpa + — [81C| - 848CaCf + 424CFTFNf\ c? + 0{o?) 



6 
+ 



15552 



7ir^(a) 
71^^ (a) 
(a) 



2 144 
2P^^ + 144 



103C| - Z%8CaCf + 104CFr^^A^/ + 0{c?) 



\cFa + ^ \31%CaCf - \\2C\ - \28CfTfN^^ 0? 



+ 



243 



27 

(5184C(3) + 20920) C\Cf - (15552C(3) + 8528) CaC\ 
- (10368C(3) + 6256) CaCfTfNj + (10368C(3) - 560) C|. 
+ (10368C(3) - 6824) CITfNj - SmCFT^Nf] + O(a^) 



72^^^ (a) 



+ 



4^ 1 
3^-" + 27 
1 



56Cl 



188CaCf + MCFTFNf\ 



243 



(7776C(3) + 4264) CaC\ - (2592C(3) + 10460) C\Cf 
+ (5184C(3) + 3128) CaCfTfNj - (5184C(3) - 280) C\ 
- (5184C(3) - 3412) C\TfN^ + 448CFr|Ar/] c? + O(a^) 



(3.15) 



where those for 7j^^(a) and 7|f^(a) are only given to two loops. Clearly the diagonal anomalous 
dimensions of each sector, including those with total derivatives, are the same as the correspond- 
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ing non-total derivative operator. Equally the mixing of the second row of W3 is trivially related 
to that of the first row of W2- This is a reassuring observation. As further checks on the results, 
since we have used the algorithm of [36] to determine the renormalization constants, therefore, 
the double and triple poles in e are predetermined by the renormalization group equations of 
(|3.10p and (|3.1ip . We note that our results are consistent with those constraints. 

We now focus on the renormalization structure of the operator correlation functions. These 
are Green's functions of operators rather than of fields but like Green's functions of fields they 
have an associated renormalization constant after the constituent bare operators have been 
replaced by their renormalized versions, taking into account any mixing. We denote these addi- 
tional renormalization constants by Z^*^^ and follow the quark current correlator renormalization 
formalism of [26 ^ 127 ^ [28]. This renormalization constant appears as a contact term rather than 
as a canonical multiplicative renormalization constant that one normally expects in the renor- 
malization of a Green's function involving only fields. In the case of operator correlators where 
there is no mixing of the constituent operators, the relation between bare and renormalized 
correlators has been given in, for example, ^26 1 \27 \ [28]. For completeness, we give the form for 
the tensor current correlation function. It is 

n5f(,) = + ;.-(z-)^n-g(,) (3.16) 

where we have included the subscript label deriving from the Lorentz tensor decomposition since 
there will in principle be a divergence for each projection. This relation is the basic form used 
in our computer algebra setup and the automatic Feynman diagram renormalization procedure 
of [36] is easy to extend and encode in Form for this case. Given ()3.16p it is straightforward 
to derive the renormalization group equation satisfied by the renormalized correlation function. 
Applying (|3.9p to (|3.16p we have 

= A^l^njf (g) + 27^(a)nj;f (g) - q%f (a) (3.17) 

where the correlation function anomalous dimension is formally given by 

7(4 («) = 

Although part of our ultimate aim is to provide the finite parts of the amplitudes we will also 
determine these anomalous dimensions to O(a^) inclusive. Indeed we use the results, such as 
()3.18p as consistency checks on the explicit renormalization constants. Aside from being gauge 
independent expressions, the double and triple poles in e are again determined by the lower 
order simple poles. Moreover, if one had not taken the issue of mixing into account this internal 
consistency check would in fact fail. 

Given this form for the quark current correlators, S and V, |26[ [271 128] . and now also T, 
their extension to operator correlators where there is mixing is subtle. Therefore, we highlight 
the structure for two cases which are {V^jl^a} and {^3,1^3}. The resulting renormalization 
group functions for the remaining cases, {V, W2}, {W2, W2} and {T, T2} can be readily deduced 
from the final relations by appropriate relabelling and ignoring irrelevant terms which, say, do 
not occur for the W2 sector. For each case we consider, the key is to simply write down the 
forms analogous to p.l6p but including all possible consistent mixings and all possible operator 
combinations in the correlation functions for that sector and ensure the equation is dimensionally 
consistent. The reason for this is that the relation between bare and renormalized correlation 
functions are entwined in an intricate way. As the {V, W3} case is simple since only one of the 
constituent operator undergoes mixing, we illustrate this by giving the three renormalization 



d 
da 



Z. 



T,T 
(i) 



(3.18) 



13 



definitions explicitly as, 



nr.; 



n 



n 



^11 ^^o(i) 



yV,b 



^22'^o (i) 



(3.19) 



where the mixing matrix elements of ()3.7p appear and the factor of (f' multiplying Z^^^-^ derives 
from the dimensionality of the actual correlator in question to ensure a dimensionless renormal- 
ization constant. The next stage is to apply ()3.9p to each of equation and then rewrite the full 
set without any bare correlators. This is algebraically tedious but it is best to start with the 
final equation since it is similar to (j3.16p whence 



= /^^<f-3(^) + (7^(a)+73t(«))n;:f-3(,) 



/ 2n2 V,ddWzi X 

W ) 7(4 (a) 



(3.20) 



with 



(a) 



yV,ddW3 



(3.21) 



Considering the second equation of (j3.19p next, after the application of ()3.9p it is necessary 
to rewrite both bare correlators which now occur back in terms of their renormalized versions. 
Throughout this and all our other similar manipulations, we always associate terms with powers 
of the momentum as contributing to the correlator anomalous dimension. Hence, we have 



= ^^<r('^) + (7^(«)+72^^(a))<f^3(,) 



+ 723 («)n(/) '{q) 



2\2V,dW-i, 



i) 



(3.22) 



and 



(«) 



e + /3(a)^ + 7^(a) + ^^^'i^) 



Z, 



V,dW'j, 
(0 



+ lV3'{a)Z^ 



V,ddW'j, 
(i) 



(3.23) 



Finally, both of the last equations of (j3.19p are required to complete the set of three renormal- 
ization group functions for the {V^, W3} case. We have 







d_ 

dfj, 

+ 7i3 («)nf4 '(g) - (q ) 7m («) 



(^) + (7^(a) + 7i?(a))<f3(^) + ^^^^l 



V,dW3 



(3.24) 



where 



7(4 («) 



- e + /3(a)^ + j'^ia) + 7i?(a) 
+ 7l^na)4r' + 7ll'(a)4f'^'. 



^(0 



(3.25) 



For the {W3, W3} case the derivation follows parallel lines though the starting point is a more 
complicated set of six relations between bare and renormalized correlators as there is mixing in 
both inserted operators. These are 



(q) 



\q' 



2t 
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'11 ^13 ^^o(i) 



dW3,dW:i 
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2e 



^11 ^22 ^^o(j) 



(Q) + ^12 ^22 ^^o(i) 



(9) 



, 7VI/3 7Vl/3TT9vy3,9^vy3/„^ , 7W3 7^/31-^^3,99^3 /„n 



^3 ^22 ^^o(i) 
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12 ^23 ^^o(i) 



^23 ^^o(j) 
W3 7H'3tt9<9W-3,9<9W3 



(,) + zr^'z^im- 



2;W3,ddW3 I ^2\i 



2e 



yW3 yW3rfW3,ddW3 

^11 ^33 ^^o{i) 



dW3,dW3 



(q) 
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+ ^13 ^33 ^^o(i) 

^9Vy3,aP^3(^2)3 
2 



+ M 



2e 



n^^3,a^3(^) + 2<3^^^3n^vv'3,aaH^3(^) 
^ (^2Tf nf^J^3.aaH^3(^) 



n' 



(9) 



^c»Vy3,99H/3(^^2^3 



2e 



7VK3 7H/3TT9W^3,9aW3/ ^ , ^W3 ^W3Tjd9W3 ,ddW3 ( ^ 
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991^3, 9<9W"3 
o{i) 



(i) 



(3.26) 



It is worth noting that our choice of the upper triangular form for the mixing matrix in fact leads 
to a simpler renormalization group equation derivation from the point of view of disentangling 
the relations to produce equations without bare correlators. Again for this sector it is best to 
begin deriving the full renormalization group equations from the final equation of (j3.26p and 
then systematically move to the row immediately above in the matrix. As this exercise is equally 
as straightforward though more tedious than the {V^, W3} case, we merely record that the final 
renormalization group equations are 







,2^3^,W3,iy3/ 



+ 72^3(„)n^H^3,aai^3(^) + 7!rn«)n^r^'"^^3(g) - (g^)^(^7''^''^^(a) 
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^ d -^ddW3,ddW3 



(q) + 273^^ (a)n 



TddW3,ddW3 



where the correlator anomalous dimensions are 
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(3.27) 



7(i) (a) 



6 + /?(a)^ + 27iT^(a) 



^^"3,^/3 



+ 271^3(^)7^^3,91^3 + 271^3 (a)Z.^3,99iy3 



7W3,9VK3 



e + /3(a) ^ + 7i?(a) + lV2'ici 
+ 7l^^(a)Z(7^"^^^ + 72^^(a)^('j;^'"'^^ + 7llna)4) 



914/^3, 9914/^3 



z 



14'3,99Vy3 
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(3.28) 



Comparing these final forms with the original relationships (|3.26|) an evident pattern emerges in 
the final renormalization group equations. Not all the original bare operators of p.26p appear 
in the corresponding equation. However, this is partly because the transformation to (|3.27|) 
involves the non-trivial entanglement alluded to earlier but in such a way that no information 
is lost. In practice there are cancellations in the derivation in such a way that the coefficient 
of certain off-diagonal elements is zero. Indeed given the upper triangular form of the mixing 
matrix and the final forms ()3.27p . one could have been tempted merely to write these down 
without derivation. 



For completeness, we close this section by recording the formal definitions of the remaining 
operator correlation function anomalous dimensions. As is apparent from comparing with their 
{y, M/3} and {VF3,VF3} counterparts there is a consistent correspondence between the terms of 
the anomalous dimensions and the form of the renormalization group function itself that means 
we only record the anomalous dimensions themselves for brevity and as an aid to checking the 
renormalization group equations. We have for those cases involving W2 
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and 
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Finally, for T2 we have 
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(3.31) 



4 Results. 



We now turn to the mundane task of recording all our results for the operator correlation func- 
tions. These are broken into subsections where the first named operator of the title corresponds 
to the operator of Figure 1. In each section, we provide the finite renormalized amplitudes 
with respect to the various projections and then the associated correlator anomalous dimension. 
In recording the finite parts of all our correlators we show the overall dimension of the amplitude 
by explicitly factorizing off the overall dependence which is not the same for each sector. We 
note that we included the powers of q'^ in the contact term of the relation between bare and 
renormalized amplitudes in the basic relation in order to identify the operator correlator anoma- 
lous dimensions in analysing the renormalization group structure. For certain sectors due to the 
nature of the total derivative operators the explicit form of some amplitudes appear in an earlier 
subsection since there is a clear relation to 0{a^) with the explicit value of the amplitudes. In 
certain cases, such as the vector correlator, the vanishing of the correlator projection to three 
loops is actually an all orders feature due to symmetry. In other cases where relations hold 
to 0{a^) this may be valid to all orders but we make no assertion beyond the order we have 
calculated to. Finally, for completeness as well as for comparing conventions, we also display 
results for {S, S} and {V, V} which are in agreement with, [261 [271 EH] . With d{R) the dimension 
of the quark representation and 




we have: 



4.1 Scalar-Scalar. 
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4.2 Vector- Vector. 

<f (9) = <lXi^{a) (4.5) 
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4.3 Tensor- Tensor. 
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4 Vector- Wilson 2. 
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5 Wilson 2- Wilson 2. 



njj;^'^^(g) 



(,^)^nj;-^^(a) 



(4.14) 



nj^^^'^^(a) 



d(ii;) 



12 



5 '''-^^ 13500 



17533 473 



8 



+ £ + —i 

225 15 



+ Cf 



/419327 1816 69266, 
+ ^C(3) + rr^i 



V 303750 
16 



135 



10125 



5 ^ 675 135 



+ 



+ 



16 



C(4) 



15838 



12 



C(3)^ 



675 
8963 



/ 12235087 
V 455625 



1350 

2606 



C(3)-4C(5) 



135 



C(3) 



135 

32 



TFNf 

3541817 
~ 1215000 

£^]Ca 



399953^ 
20250 ' 



C(4) + 24C(5) + 
5 ^ ' ^ ^ 20250 



1448 128 
225 135 ' 



+ 0(a3) 



(4.15) 



n 



W2,W2 
(2) 



(a) = d{R) 



225 15 



Cf 



8,,„, 15683 986, 16,0 

-C(3) H 1 £^ £^ 

5^^ ^ 20250 675 45 



/ 1685066 3632 ,, , 48404, 

1 C(3) H £ 

V 1366875 405 ^ 10125 



19 



S2 

+ —C(3)£ + 
15^^ ' 675 



1084 ,2 V ^ 



+ 



+ 



64 
405 



15^ ^ 



29276 _ 8 _ 7819793 
2025 '~ b'^^ '~ 2733750 



46547^ 
3375 ^ 



- 5«^>'- 675 

/ 8792588 5788 



3061^2 _:^^3)^^ 



V 455625 



405 



C(3) 



176 

405 

64 



+ 



9008 
2025 



■r + 



256 
405^ 



44597 

Cf] a^] + O(a^) 



(4.16) 



d{R) 



17 2 

+ —i - Cf 

225 15 



34439 



/464 

+ ic(3). 

/ 918157 
V 15000 " 



6750 

6263527 128096 



638^ _ ±^2 
225 15 



303750 
1916 



10125 



+ 



+ 



675 
17492 



675 ^ 5 ^ 3^^ ^ 10125 



15 ' 675 135 J 

, 32^,^, 436,, , 28297949 
16C(5) + -C(4)--«3)-^^ 



160709 
10125 ^ 



128J 



12^ 
135 



t'\CF 



+ 0(a3 



(4.17) 



n;yf^^'''^(a) 



d{R) 



55 



- 8C(3) + 2e 



304 



C(3) 



3701 44 32 , , 4 , 



+ 



/ 44215 



V 324 
+ |80C(5) 



^<(3)-fc(5)+4U 



11 



3-C(3)^ + f^')CA 



148 



C(3) 



143 



Cf 



+ 0{a^ 



(4.18) 



nf,^^'^"'^(a) 



n 



^1) 

aw2,dW2 



(3) 



(a) 



0(a^ 



2n 



14^2, a Wa 
(3) 



(a) + 0(a3 



(4.19) 



d{R) 



I + 4^[(259200C(3)-65603)C^ 

+ (325498 - 518400C(3))Cf + 22612rFArj] a^l + O(a^) 



20 



d{R) 



d{R) 



+ (345738 - 518400C(3))Cf - 7828rFAr/] a^l + 0{a^) 
+ (259200C(3) - 184754)Cf - 26276TFAry] a^l + O(a^) 



+ 0(a3 



















2 

- + 2CFa + 

O 


Cf 
18 


[133Ca 


- 18Cf 


- A'iTFNf] 




= 7(2) (a) = 


Cf 
9 








7^3^'^^^ (a) 




4 

- + ACpa + 


[133Ca 


- 18Cf 


- \\TfN^\ c? 




+ 0(a3) . 











(4.20) 



4.6 Vector- Wilson 3. 



(0 



(4.21) 



n(i) '(a) 



n(2) '(«) 



= d(i?) 



.675 

+ Cf 



■^-1C(3) + ^^' 
6480 15^^ 27 



(ic<3) 



4070273 7267 



2624400 14580 



' + ^C(3)^-— 
45^^ ^ 243 



n 



(2) 



= d{R) 



+ 



+ 



a = 



— + — ^ + 
,27 27 

Y152 , , 

/ 44437 
V 5832 



/ 48524449 1364 
V 10497600 ~ los" 
4051 



C(3) - -C(5) 



+ 



2916^ 45 
„ _ 1118783 
^ 2624400 



972 
^C(3,- 



187 
19440 



i\CF 



\(^) — " 2 (2) 



2^(1) 



+ Cf 



18 



19 4 , , 1 

C 3) + -I 

36 9^ ^ 9 



3719 22 „ 16 _ 2 _ 

27^ j T.iV, 



-^f + -^C(3)^ 
1458 27 27^^ ^ 



+ 



fc(3,-|c(.) 



+ 



+ 0(a3) 



81 
44 
27 

9 27^^ ^ 



27 

Sc(3K.ii.^)c. 

— ^) Cf 
18 ' 



(4.22) 



21 



7^f (a) = - 75r(«) = ^(^) 



V,W3, 



1 13 ^ 
45 + 162^-" 



7(1) (a) 



+ 



7(^'r(a) 
1 



d{R) 



1 



2^(1) 

1 „ Cf 



116640 
(a) = 



[27062Ca + 239Cf - ^imTpNf] 
1 



v,ddW3 
2^(2) 



(a) 



27 + 9GFa + Y^[37CA-6Ci.-12TFiV;]G2 



+ 0(a2 



+ 0(a3) 



(4.23) 



4.7 Wilson 3- Wilson 3. 
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4.8 Tensor- Trans versity 2. 
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Finally, we note that in addition to the various checks we have mentioned so far, our Form 
code was written in such a way that only the Feynman rules for the operators and projectors 
needed to be input. The Mincer integration code and its interface with the Qgraf set of 
Feynman diagrams forms the same central block module of the programme. In this way our 
approach was designed in order to minimize the potential places where errors could creep into 
the overall computer algebra computation. In this respect we have not in fact derived new 
Feynman rules for the parent operators V , T, W2, W3 or T2 but imported those used in the 
progammes which underlay the results of [21^ I22j. For the remaining total derivative operators 
it is evident from the consistency, say, in the relations of their anomalous dimensions with those 
without the derivatives that their Form Feynman rule module is not inconsistent. 



5 Tensor current i?-ratio. 



Our final exercise is to derive the i?-ratio for the tensor current to complete the evaluation for all 
the quark bilinear currents. As for S and V it can simply be derived from the current correlator 
by using 

where e indicates the usual shift away from the real axis to avoid ambiguity. Unlike S there 
are two Lorentz tensor components and for the moment we assume there are two i?-ratios. The 
case V has in principle two similar channels but due to gauge symmetry there is no contribution 
in the longitudinal piece of the decomposition. With this definition and our results (14. 9|] and 
(ITOD we find 
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where 




(5.3) 



From these it is evident to see that there is a simple relationship to three loops between both 
channels which is 

(5.4) 



^f2)(^) 



2i?fi)(.) + O(a^) . 



Whilst the full expressions for ()4.9p and ()4.10p are different and do not satisfy an analogous 
relation, the behaviour of the I terms do which is the origin for the result ()5.4p . Consequently, if 
we now include the Lorentz tensors of the projection basis we can write down the Lorentz tensor 
dependence of the i?-ratio as one would have derived it directly from ^J^^i^vw^^'^^ 
had the problem of focusing on scalar amplitudes in order to perform the Mincer calculations. 
Therefore, we have 
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where we have introduced the common tensor structure 
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The appearance of this tensor structure is akin to that for case V where the longitudinal piece 
is absent. Put another way this form would have emerged directly if we had chosen our Lorentz 
tensor basis in a more erudite fashion. For completeness, we have numerically evaluated the 
amplitude for the colour group SU{3) similar to Appendix B. We have 
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For example, to see the convergence behaviour in relation to the expressions for S and V for 
three quark flavours when s = fi'^ we have 
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0.333333 + 1.037037a + 1.25359402] + 0{a^) 
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6 Discussion. 



We conclude with brief remarks since the main goal of the exercise to determine the finite 
parts of various operator correlation function to 0{a'^) in the MS scheme has clearly been 
achieved. It extends the work of [261 113 [2H]- One novel feature was the need to properly 
account for the operator mixing into total derivative operators for the flavour non-singlet twist- 
2 operators used in deep inelastic scattering. The mixing matrix has been deduced for several 
low moments but to extend these to moments n > 4 for arbitrary n to even two loops would seem 
to be excluded at this stage. For instance, the calculational machinery on a par with Mincer 
is unfortunately not available. Whilst the main obstacle is the inability to disentangle the 
relations between counterterms one way through could be to embed the operators in higher leg 
Green's functions. Whilst this, in principle, will give more relations between the counterterms 
there is again the problem of lack of calculational machinery. Indeed with more legs with 
independent momenta any nullification of external momenta has the additional potential problem 
of introducing spurious infrared singularities. These would have to be properly treated using, 
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say, infrared rearrangement to be confident in the correctness of the final counterterm relations. 
However, since the main problem here was motivated by the need to provide only low moment 
flavour non-singlet information for lattice computations, this is a problem which is left for future 
consideration. 
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A Projectors. 

In this appendix we record the explicit forms of the tensors into which the various correlation 
functions are decomposed. For each sector wc also record the matrix A4^^ used to project out 
each individual component of the decomposition. The matrix Ai^j^i is derived by first constructing 
the matrix A/^j where k and I label the projectors, which is defined by 

^kl - '^{k){n^...,Jir,^\vi..M^.}\^)l^{l) W l^-tj 

where there is no sum over the i and j. The elements of this matrix are polynomials in the 

dimension d due to the contraction of the Lorcntz indices. Finally, M"^^ is the inverse of A/"*^. 
Once Al*-' is specified then to project out, say, the fcth piece of the tensor correlation function, 
one multiplies it by the projector 

1=1 

where there is no sum over the labels {ij}. The method we have used to construct the tensor 
basis, which of course is not unique, is to first write down the complete set of tensors built from 
the metric, ry^j/, and the momentum, q^, which have the same number of free indices as the 
operator correlation function of interest. Each of these independent tensors is then multiplied 
by a different label and then the Lorentz symmetry properties of the two operators in the 
correlation function are enforced on the sum of all independent tensors. This provides a set 
of linear equations for the labels which is fewer in number than the total number of original 
labels. Solving these equations reduces the number of independent labels, and hence independent 
combinations of the individual tensors, producing the tensor basis as enumerated in Table 1. 
Therefore, it remains to list the relevant explicit expressions for the various sectors as: 



A.l Vector- Vector. 



/^{innMW - ^2 ' ^"{2){^,\v}Vl) - ^2 y^-"^) 

M^'^ = i ° ^ (A 4) 

{d - I d-1 ^ ' 



A. 2 Tensor- Tensor. 
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T T 



/ \ QuQp iiyia ifiip 



Qp.qp 



M 



4(d- l)(d-2)g2 1-2 d 



3 Vector- Wilson 2. 



1 ..{^^xyQj 



:,v,W2 

{2){n\ap} 



(?) 



,qpqpQa_ 



M 



V,W2 _ 



d 



2d{d - l)g2 10 2 



4 Vector- Wilson 3. 



yV,W3 
{l){p.\apX} 



{2){^i\crpX} 



(?) 



- {d + 2) [rif,aqpqx + rji^pQaqx + V^xqpqa] ^ 

+ 2(d + 2)M^ 

[vapqij.qx + v<7xqp.qp + Vpxqp.qa]-^ - {d + 2) ^^^^^y^ 



M 



V,W3 



1 - 1 



3(d-l)(d-l)(d + l)(g2)2 I - 1 3ci + 4 



5 Wilson 2- Wilson 2. 



(l){/il/|o-p} 

3W2,W^2 

{2){pu\ap} 

(3){|Ui/|o-p} 



(?) = f]txar]vp + r/^pryi.a - -^Vixi^Vap 

(?) = - ^'^iJ.uVap + [Vp-vqaqp + Tlapq^iqy] ^ 



(?) = [vnaq^qp + ri,_,pq„qa + v^aq^qp + v^pq/xqa 



- 4 



qpq^qaqp 



q" 



M 



W2,W2 _ 



4(d-l)(d+l)((i-2)(g- 



.2A2 



/ 2(d - 1) 4 - 2(d - 1) \ 

4 4d - 4 

^ - 2(d - 1) - 4 (^2 + d - 4) y 
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A.6 Wilson 3- Wilson 3. 



yW3,W3 
(l){/iI/CT|pAV)} 



■yWs^Ws 

{2){fj.i^a\p\^} 



r, 



W3,W3 

(3){/ii/CT|pAV>} 



flliuVcrpriXi' + VfiiyllaWfrip + llt^uVaijIlpX + rjtia'nup'nXtp + VnaViyXVpip 
+ Vficrllu^llpX + VfipVuallXtlj + IJuXVuallpilj + V^^-q^a'TlpX 

{d + 2) . 

^2 — mvrjapqxQip + ViJ,uV<7xqpqi> + "n^v-na^qxqp 

+ Vnariupqxqi^ + VnaVuxqpqi^ + VnaVi^^qxqp 

+ Vp.pVuaqxqi> + Vp.pVxipquqa + VfixVuaqpq^ 

+ Vt^xVpipQiyQcT + Vti'PVvaqpqx + Vti^Vpxq^qa 

+ ijuptix^q^qa + Vvxvpipqfiqcj + Vutpiipxqpqa 

+ riapr]x^q^iqy + VaXVp^qi^qu + Va^jVpxqi^qu] 



+ 



2{d + 2) 



(g2)2 



bip-uqaqpqxqf + ViJ-aq^qpqxq^ + Vvaqp.qpqxq^ 

+ -npxq^q^q^q^ + ■qp^qp.q^q^qx + rix^q^qi^qaqp] 



(d + 2)2 ^ 

H — r^TT^ [nppququqxqip + Vtixqvqaqpqi, + Vp.'^pqvqaqpqx 



8{d + 2) 



+ rii,pq^q„qxq^ + Vvxqiiq^qpq^ + Vui^qfiqaqpqx 

+ Vapqp,quqxqi^ + Vaxqp.quqpqi, + r]„^qpquqpqx\ 



(g2p — qiJ.q'^qTqpqxq-ii) 

riixpriuxriailj + rif^pTiu^^riaX + »7^a»7z^p'7ctV 

+ Vnxriwipriap + 'nn^r]j,pr]ax + iJi^^rji^xriap 
2 ^ 

- [iip.ur](rpqxqtp + 'n^u'naxqpq'4, + 'niJ.u'nai(,qxqp 

+ Vfj.aVupqxq'il, + 'n^a'nvxqpq^, + 'n^^'q^^qxqp 

+ VppVuaqxqi^ + VppVxi^quqcr + VpxVi^aqpqi; 

+ r]^,xr]fni,qvqa + r]i^i>riuaqpqx + rj^^Vpxquqa 

+ VupVxrpqp.qa + VuxVpipqiJ.q<j + rjv^Vpxqixqa 
+ VapVxfq^qi^ + rjaxVp^qp-qv + VaipVpxq^qu] 



+ 



+ 



(q 



2^2 



[vpuqaqpqxqi^ + Vpaquqpqxq^, + 'n^aqfiqpqxqi, 



2{d + 2) 



(g2)2 



Vpxqpqyqaqi, + rip^q^q^qaqx + riX'^qpqvqaqp] 
bli^pquqaqxq^ + Vnxqi^qaqpqip + Vnipqi^qaqpqx 



+ r]vpqp,qaqxqi> + Vuxq/^qaqpq^ + Vu^q^tqaqpqx 

+ Vapq^quqxq^ + Vaxq^quqpq^ + Vai^q^quqpqx] 

16(d + 2) 

(g2-)3 QnQ'^Q'rqpqxqi' 

^ [vp-uVpxqcrq^p + Vp-uVpipqaqx + Vixi^vx^qaqp 

+ Vfj-aVpxquq^ + Vp-aVptpqi^qx + ViJ.am^qvqp 
+ VufjVpxq/xqf + VuaVp^pqp.qx + v^aVx^q^qp] 

{d + 2) 



{q 



2\2 



[viJ.i^qaqpqxq^ + Vp.aqi^qpqxq^ + Vi^^qt^qpqxq^ 
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+ rjpXQfiQuqaq^ + Vpi^qixQuqaqx + rixi^q^quqaqp] 

+ -j^^qij.qi^qaqpqxqf 

^mZUpm^^^ = ^ [^/w^i^A^ag^ + Vi^pVu^q^qx + Vt^pVaxq^q^ 

+ Vf^pVai,quqx + Vt^xViypqcrq-ip + Vt^xViyH^qp 

+ Vt^XVcrpqiyqi, + Vt^XVaipq,yqp + Vt^i>VupqaqX 

+ riti^ri^xquqp + ri^^rj„pq^qx + rj^^rj^xqvqp 

+ r]yp'naxqp.qi, + Vi^pria^q^qx + v^xVapq^qi^ 
+ VuxVaipqfiqp + Vu^VcrpqfMqx + Vi^iPVaxqtiqp] 

2 . 

[Vpvqaqpqxq-4, + Vpaquqpqxq^p + Viyc7qfj.qpqxqt 



(g2)2 

+ Vpxqnqi^qaqi> + Vp^pqi^q^qaqx + m^q^qi^qaqp] 

4 ^ 

- -^^2)2 uip.pqi^qTqxq^ + "nixxqvqaqpq^ + Vp.Huqc7qpqx 

+ rj^pQ^Qaqxq-^ + Vi^xqixqaqpqip + r]„^qfiqaqpqx 

+ r]apqp.quqxqil, + Vaxqfiquqpq-,p + Vatjjququqpqx] 
2(d+14) 

^ — (-g2-)3 (it^it^iTqpqxq^ (a.is) 



W3,W3 



1 



18(^2 - - 2)(d + 2)2(d + 3)(g2)3 



/ 2(7d + 18) 
-6(d + 2)2 
- 2{7d + 18) 
\ 6(d + 2)2 



-6((i + 2)2 
3(ci + l)(d + 2)2 

6(d + 2)2 
-3(d + l)(d + 2)2 



- 2{7d + 18) 
6(d + 2)2 



6(d + 2)2 \ 
■3(d+l)(d + 2)2 



2(lld2 + 50d + 48) -2(d + 6)(d + 2)^ 



-2(d + 6)(d + 2)2 



d(d + 5)(d + 2)2 y 
(A.14) 



A. 7 Tensor- Transversity 2. 



{l){p.u\apX} 



iq) = 



V, 



T,T2 

{2){pu\apX} 



(q) 



:>T,T2 

{2,){lJ.v\apX} 



{q) = 



'^'{A){pu\apX}'^l) 



Vp.aViypqX + VixcrriuXqp - VlxpVi^aqx - V/xXriuaqp 

+ Vip.pqvqaqx + Vp-xqi^q^qp - Vvpqij.qaqx - v^xq^qaqp 

+ ^ri„aqnqpqx - ^rjp^q^qpqx] 

VfiaVpXqu - Vuar]pXqfi 

+ [Vupq^qaqx + Vuxqfiqaqp - ri^pq^q^qx - Vpxq^qaqp 
+ dr]uaqp,qpqx - drji^^qi^qpqx] ^ 

riixpTjaXqu + ■nixXrjapqv - VupVaXqu - rjvXrjapqp. 

+ [{d + 1) {vupq^qaqx + Vuxq^qaqp - Vi^pq^qaqx - Vnxq^qaqp) 
+ ^Vuaqp-qpqx - ^Vpxrqi^qpqx] 

"np-pnvXqa - Vi^XVupqa 
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+ [Vup%(lcrq\ + VuXQixQaQp - VlJ-pQi^QaQX - VfJ-XQi^QaQp] -2 (A. 15) 



( (f{d+\) 



\ 





2{<f + 4) 

- Ad 
A{d - 2) 



\ 

- Ad A{d - 2) 

- 2d{d - 2) 
- 2d{d-2) 2(d-l)(d2-4) J 

(A.16) 



B Expressions for SU{3). 

For completeness and for practical use, we record the explicit numerical values of the various 
amplitudes for the colour group SU{3). We take the usual values for the Casimirs, Tp = ^, 
Ca = 3 and = | as well as d{R) = 3 but leave the numbers of quarks unfixed. We only 
record those amplitudes which are non-zero. The remaining ones still satisfy the same relations 
to third order which were noted in Section 4. Thus, we have 



B.l Vector currents. 



n^'^(a) 



n(i') (a) 



4.000000 + 2.000000£ + 48.867512 + 45.333333£ + 8.000000^ 
+ [1925.894130 - 100.119646iV/ + (1650.138715 - 61.022786Ar/) 
+ (565.333333 - 19.5555557V/) f 
+ (50.666667 - 1.777778iV/) a^l + 0{a^) 

3 [- 2.222222 - 1.333333^ + [1.199436 - 5.333333^] a 

+ [6.784729iV/ - 38.534112 + (2.459635iV/ - 42.361758) i 



+ {i.mnmf - 29.333333) f 



+ 0(a3) 



= 3 



- 0.444444 - 0.666667^ + 



0.698484 - 2.074074^ + 0.888889r 



+ [27.577316 - 1.114284A7/ + (22.743851 - 2.243433A7/) i 
+ (18.172840 - 0.592593iV/) f 



n5f(a) = 3 



+ (2.469136 - 0.197531A(^) 



2.222222 + 1.333333^ + 



3.640070 + 4.148148^ - 1.777778^^ a 
+ [- 32.988175 + 2.2724637V/ + (- 45.487702 + 4.486867A7^) i 
+ (- 36.345679 + l.lShl^hNf) 



+ (- 4.938272 + 0.395062Ar/) . 



+ 0{c?) 



(B.l) 



B.2 Wilson moment n 



nJ5'^2(a) = 3 



0.480000 + 0.200000^ + 



- 5.578236 - 2.802963^ - 0.71111W^ 



+ [- 88.800839 + 11.700261Ar/ + (- 56.861337 + 7.1251127V/) ^ 
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n;rf^^(a) 



+ (- 15.116049 + 1.566420iV/)£^ 



+ (- 0.921811 + 0.158025iV/) 
0.408889 + 0.133333i' + 



+ 0{^) 

3.597014 - 1.947654^ - 0.474074^^ 
+ [- 57.728474 + 8.008477A[/- + (- 39.368555 + AmQmNf)i 
+ (- 10.231001 + 1.070617iV/)£2 



+ (- 0.614540 + 0.105350A(/^) 



0.075555 + 0.133333^ + 



+ 0(a3) 



4.238377 + 3.780741^ + 0.711111£^ 
+ [75.027286 - 10.992767iV/ + (56.696221 - 6.724712Arj) ^ 
+ (21.758025 - 1.892346iV/) ^2 

+ (0.921811 -0.158025Arj)^^l a^l + O(a^) 

3 [1.111111 + 0.666667£ + [- 0.599718 + 2.666667^] a 

+ [19.267056 - 3.392365Ar/ + (21.180879 - 1.229818Ar/) £ 



+ (14.666667 - 0.888889iV/; 



+ 0{a^) 



(B.2) 



B.3 Wilson moment n 



n^f(a) 



n^i) (a) 



3 [0.045926 + 0.022222^ + [0.020544 + 0.098765^] a 

+ [0.917043 - 0.131672Arj + (0.838448 - 0.047350Arj) £ 
+ (0.497942 -0.030178Arj) a^l + O(a^) 

3 [0.074074 + 0.037037^ + [- 0.008626 + 0.148148^] a 
+ [1.211681 - 0.196695iV/ + (1.176715 - 0.068323iV/) £ 



+ (0.814815 - 0.049383iV/) 



+ 0(a3 



(B.3) 



nf[^.i^3(a) 



ng'^^(a) 



0.001151 + 0.000265^ + - 0.037289 - 0.020408£ - 0.003527^^ 
+ [- 0.344878 + 0.084926iV/ + (- 0.203993 + 0.046931iV/) I 
+ (- 0.050228 + 0.009964iV/) 



+ (0.000131 + 0.000784iV/) 



+ O(a^) 



0.004528 - 0.001587£ + 



0.100640 + 0.049139£ + 0.008818^ 



+ [1.022184 - 0.224027Arj + (0.551061 - 0.122691iV/) i 
+ (0.109786 - 0.024512iV/) f 



+ (- 0.000327 - 0.001960Ar^) ^3 



0.010207 + 0.003439£ + 



+ 0{a^ 



0.069547 - 0.035555^ - 0.00964ir 
+ [- 1.354970 + 0.153442Arj + (- 0.921737 + 0.095085A(y) £ 
+ (- 0.233863 + 0.020396Ar/) ^2 



+ (- 0.017201 + 0.002143iV/) 



+ 0(a3 



- 0.008806 - 0.003968£ + 



0.063950 + 0.031660£ + 0.010935r 
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n 



W3,ddW3 
(1) 



n' 



dW3,ddW3 
(1) 



(a) = 



[1.580679 - 0.12m7Nf + (1.078542 - 0.084549A(^) £ 
+ (0.309941 - 0.019979iV/) ^2 



+ (0.025932 -0.0024307V/) c 

3 [- 0.001481 - 0.000617^ + [0.002510 - 0.001641^] a 
+ [0.048778 - 0.003811iV/ + (0.009906 - 0.0018187V/) ^ 



+ (- 0.011698 + 0.000472iV/) r 



0.012840 + 0.00432W + 



+ 0(a3 



0.102610 - 0.052455^ - 0.013169£^ 
+ [- 1.653254 + 0.226264iV/ + (- 1.103477 + 0.137837iV/) £ 
+ (- 0.272497 + 0.029267Ar/) £2 



+ (- 0.017071 + 0.002926Ar/) ^3 



+ 0{a^ 



0.013333 - 0.005555£ + 



0.154951 + 0.077860£ + 0.019753^ 



(a) = 3[ 



+ [2.466690 - 0.325007Ar/ + (1.579482 - 0.197920iV/) i 
+ (0.419890 - 0.043512iV/) 
+ (0.025606 -0.004390iV/) a^l + O(a^) 

3 [- 0.002963 - 0.001235£ + [- 0.002132 - 0.005487^] a 

+ [- 0.055664 + 0m7575Nf + (- 0.046580 + 0.002631iV/) £ 
+ (- 0.027663 + 0.001677iV/) a^l + O(a^) 

3 [- 0.002401 - 0.001029^ + [0.003129 - 0.003018^] a 
+ [0.039866 - 0.001957iV/ + (0.000510 - 0.001236iV/) £ 



+ (- 0.020450 + O.OOlOOeiV/) £■ 
0.013759 + 0.004733^ + 



+ O(a^) 

0.103229 - 0.051084£ - 0.013169^^ 
+ [- 1.644342 + 0.2244107V/ + (- 1.094081 + 0.137255iV/) £ 
+ (- 0.263695 + 0.028733iV/) £'^ 
+ (- 0.017071 + 0.002926iV/) a^j + 0{a^) 

3 [- 0.004801 - 0.002058£ + [- 0.000893 - 0.008230^] a 

+ [- 0.073488 + 0.011283iV/ + (- 0.065373 + 0.003796iV/) £ 



+ (- 0.045267 + 0.002743Ar/) ^2 



(B.4) 



B.4 Trans versity moment n — 2. 



n 



T,T2 
(1) 



(a) 



njf (a) 



0.222222 + 0.333333^ + - 0.349242 + 1.037037£ - 0.444444^^ 
+ [- 13.788658 + 0.557142iV/ + (- 11.371925 + 1.121717Ar/) £ 
+ (- 9.086420 + 0.296296iV/) £'^ 



+ 0{a^ 



+ (- 1.234568 + 0.098765iV/) £^ 

0.740741 + 0.222222^ + [2.335073 + 0.691358^ - 0.296296£2 
+ [5.429324 + 0.264127iV/ + (- 7.581284 + 0.74781 liV/) £ 
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+ (- 6.057613 + 0.197531iV/) ^2 

+ (- 0.823045 + 0.065844Arj)£3l j^. o{a^) 



(B.5) 
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